Calculation of the Branching Ratio of B — > h c + K in PQCD 
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The branching ratio of B~ — > h c + is re-evaluated in the PQCD approach. In this theoretical 
framework all the phenomenological parameters in the wavefunctions and Sudakov factor are priori 
fixed by fitting other experimental data, and in the whole numerical computations we do not intro- 
duce any new parameter. Our results are consistent with the upper bounds set by the Babar and 
Belle measurements. 

PACS numbers: 13.20.He, 12.38.Bx, 13.60.Lc 



I. INTRODUCTION 



h c of 1 P\ is the spin-partner of r\ c in the charmomium family, however, it seems to h c behaves very differently 
from other family members. Searching for h Cl as well as hf, is a challenging task for both experimentalists and 



21, and later 



theorists of high energy physics h c was first seen at the CERN Intersecting Storage Rings (ISR) 
the E760 group at Fermilab reported observation of h c by studying the pp — > h c — » ir°J/ip ^|. However, the 
E835 group claimed that they did not see h c in this channel. Instead, E835 reported observation of h c in another 



channel pp 



hr 



IVc 



111 



J]. Recently, the Babar collaboration set an upper bound for the production 

By contrast, the Belle 



rate of h c as 3.4 x 1CT 6 for BR(B~ -> h c K~) x BR(h c -> J/V'tt+tt") at 90% C.L. 
collaboration reported null result on searching for h c via B ± -> hcK* 1 and gave an upper bound on the branching 
ratio as BR(B+ -> h c K+) < 3.8 x 10" 5 

ip(2S) — ► TT°h c — > 7r°7?7 c with a very small branching ratio [7|, l8(. As a comparison, so far, the BES II collaboration 
has not seen h c yet, and searching for it may be an important issue for the upgraded stage BES III. 



fjj. The CLEO Collaboration has announced observation of h c in the decay 
ratio 0,3. 



One needs to understand the experimental status about the production rate of h c . Moreover, a relatively larger 



production rate may be associated to new physics 



3,3 



ll| , therefore study on h c production may provide a chance 



to investigate effects of new physics, such as supersymmetry at lower energy scales. Of course, before invoking for 
new physics, one needs to more accurately estimate the production rate of h c in the standard model (SM). 

Suzuki suggested to look for h c at B + — > h c K + — > r )rj c K + — > j(KKir)K + and if approximately BR(B + — > h c K + ) as 



BR(B + — > XcoA" + ), he estimated the cascade branching ratio as 2 x 10 



12j. Gu considered another decay chain 



and estimated the branching ratio BR(B + — ► h c K + — > j?] c K + — > 7(A"gA' + 7r +c.c.)A" + — > 7(7r + 7r A" + 7r +c.c.)A + ) 



as 3.5 x 10" 6 



13] 



In the decay of B-mcson into charmonium is the so-called internal emission process where due to the color matching 
the process is suppressed compared to the external emission. Moreover, the non-factorizable effects would further 



change the contribution of the internal process besides a color factor of 1/3 



14[ . Therefore it seems that the smallness 



of B — ■> h c + K is natural. However, one could not reproduce experimental value of the B — > XcoK decay rate 15 1 



in the QCD improved factorization, so that the authors of 
decay. They claimed that 



15j suggested to consider re-scattering effects in B meson 



15| a larger branching ratio for B~ — > \coK~ which is consistent with data, was obtained. 
Motivated by the same idea, they applied the same scenario to investigate the re-scattering effects for B~ — > h c K~ 
case, i.e. supposing that the decay B~ — > h c K~ occurs via the the subsequent re-scattering effect of Z?i*' — D^°, 



the products of decay B~ -> D^D^ [l6|. They obtained a branching ratio (BR(B~ -> h c K~) = (2 ~ 10) x 1CT 4 ) 
which is much larger than the upper bound set by the Belle collaboration. It is very possible that the branching ratio 
of B~ — > h c K~ was overestimated in their work due to uncontrollable theoretical uncertainties, such as some input 
parameters and basic assumptions adopted in their calculation, which were comprehensively discussed in their paper. 

Even though we may suppose that we have full knowledge on the weak and strong interactions at the quark level 
and can derive the quark-transition amplitude, the most difficult part is evaluation of the hadronic matrix elements of 
the exclusive processes. In fact, at present, a complete calculation of i? ± — > h c K^ based on an underlying theoretical 
framework is absent. The perturbative QCD (PQCD) approach is believed to be successful for estimating transition 



171 ]. even though there is still dispute about its applicability 



18(. The authors 



rates of B and D into light mesons 
of ref. 19( applied the PQCD approach to study B~ — > XcoK~ and also obtained results which satisfactory comply 
with data. Therefore we have reason to believe that it is appropriate to analyze B^ — > hcK^ in this framework. In 
this work, we calculate decay rate of B~ — > h c K~ in the PQCD. 

Generally, for two-body non-lcptonic decays of B meson, both factorizable and non-factorizablc diagrams contribute 
to the _ -ever, * B" - MT", «— utions tan t _le dia gmls disappea, since 

the conservation of G parity leads to (0 [03^75 c|fo c ) = [20(. Therefore, the decay rate of B^ — * h c K^ is much 



different from that in i? ± — > J/tpK^ 



19 



2l| where both factorizable and non-factorizable diagrams contribute. 
To be more precise than the qualitative understanding, one needs to calculate the non-factorizable contribution 
where the QCD effects are accounted. Moreover, since the process is not factorizable, the convolution integral would 
involve the initial B-meson and all the two produced mesons altogether. 

Our numerical results indeed indicate that the order of magnitude of B~ — > h c K~ should be of order of 3.6 x 10 -5 , 
which is smaller than the upper bound set by the Babar collaboration 5] and slightly below the upper bound given 
by the Belle collaboration 6j. It is noticed that our result about B(B~ — > h c K~) is smaller than that estimated by 



the authors of 



12 



13 



15j and more consistent with data. 
Indeed, more decisive conclusion should be made as more accurate data are accumulated by BES III, CLEO, Babar, 
Belle and even the LHCb. 



The structure of this paper is organized as follows. After this introduction, we formulate the decay amplitude of 
B~ — > h c K~ in the PQCD approach. Then we present our numerical results along with all the input parameters in 
Sect. III. The last section is devoted to our conclusion and discussion. Some tedious expressions are collected in the 
Appendix. 

II. FORMULATION 



The quark-diagrams which contribute to the transition amplitude of B~ — > { l Pi)h c + K~ arc displayed in Fig. [T] 
As has been discussed in the introduction, the factorized diagrams Fig. [1] (a),(b) do not contribute to the amplitude 
because the the G parities of axial meson h c and axial current are mismatched and in the flavor SU(3) symmetry 



limit, the corresponding hadronic matrix element is forbidden 



Q.Q. 



(a) 





(d) 



FIG. 1: Feynman diagrams correspond to the calculation of hard amplitudes in B — > h c K. 



The effective Hamiltonian relevant to B — ► h c K decay in the SM is written as 22 1 

10 



H 



eff 



Gp 
V2 



V cb V c * s (C 1 ( l i)0 1 + C 2 (ti)0 2 ) - V tb V t * s J2 C ^)°i 



i=3 



where Ci(/x) are the Wilson coefficients and Oi are the relevant operators defined as 

0\ = (SaCp)v-A(cpba)v-A, 2 = (s a C a ) V ^ A (c/jbp) V ^A, 
^3(5) = {Sgbg)v-A y^^Plp)v-A(V+A), ^4(6) = {Sgbp) V-A /~](QpQa)v-A(V+A) ) 



O 



7(9) 



{SabajV-A^eqiqpqf}) 



V+A(V~A), ^8(10) 



\{s a bp)v -A^2e q (qf}q a ) 



(1) 



V+A(V-A): 



with a, (3 being the color indices. The explicit expressions of the Wilson coefficients appearing in the above equations 



can be found in Ref. 



22] 



We define, in the rest frame of the B meson, p, p' and q to be the four-moneta of B, K and /i c , fci(2), ^(2) anc ^ ?i( 2 ) 
to be the momenta of the valence quarks inside B(6(q)), K (s(q) ), and h c (c(c) ), respectively. Then we parameterize 
the light cone momenta with all the light quarks and mesons being treated as masslcss 



p= ^(1,1,0 t ) = (p+,p-,0t), p' = ^J(0,l-r 2 ,O T ), q = ^(1, r 2 , T ) = (q+ , q~ , T ), 
fei = (xip + ,p~,ki T ), &2 = {xip + ,0, -ki T ), k[ = (a;ip' + , x[p'~ , k' 1T ), 
fc 2 = 0^il> /+ )a ; 'iP /_ ,-k / 1T ), 9i = {yq + ,yq~, c tr), q 2 = (yq + ,yq~,-qr), 

where the mass ratio r is set as r = mh c /mB- Xi and x\ are the fractions of the longitudinal momenta of the valence 
quarks. The superscripts "+" (plus) and "— " (minus) mean that the three-momentum is parallel or anti-parallel to 
the positive z direction which is defined as the direction of the three- momentum of the produced h c . kix, k' 1T and 
qx are the transverse momenta of the valence quarks inside B, K and h c respectively. 

Therefore, we concentrate on the calculations of the non-factorizable diagrams Fig. [T] (c),(d). Then we divide the 
operators appearing in the effective Hamiltonian into two categories according to their chirality. i.e. (V — A) ® (V — A) 
and (V - A) ® (V + A). 

For the type of (V — A) g) (V =p A), the hard kernels of Fig. [1] (c) arc respectively written as 



m B 



m B 



H 



H 



a/3pa'/3'p>(P'P ' l) 



[7 M (l-75)]a'p'[*fl«7 



■(1-75) 



[(1 + 7S)]aV [idsYhp' TT _u ]2 - 
pa \ 2 2 / 



where the factor '-2' in eq. ([3]) comes from the Fierz transformation on the (V — A) <£> (V + A) operators. 
Similarly, for Fig. [T] (d), we have 



(2) 
(3) 



ff i£ W (P.P'9) = [7m( 1 -75)] pq 



igsYii - 75) 



-(fo- fo+ h' 2 ) - m 



■1v 



H 



tf3pl>0> P >(Pi'P'i<i) = - 2 [7/ 1 (l -75)U 



^7^(1 +75)- 



"7* 



(4) 
(5) 



-W2- ^2+ A2) - ™ c 

In the above expressions, indices 1,2 denote the contributions from {V — A)®(V — A) and (V — A)® (V^ + A) operators 
respectively. 

Finally, for each type of operators we obtain the decay amplitude which is a sum of the two non-factorizable 
diagrams. For saving the space, we collect them in Appendix A. The wave functions relevant to calculation is given 
in Appendix B. 

The decay width of B~ — > h c K~ is written as 



r 



G 2 F |P/| 
2 8ttto 2 r 



(6) 



with the total amplitude is the sum of Mi^(c, d) as 

M = M[ c) + M 2 (c) + M[ d) + M^. 

Pf denotes the three-momentum of the produced meson in the center-of-mass frame of B meson. All the explicit 
expressions of Mj are collected in Appendix A in order to shorten the text and focus on the physics contents. 

III. NUMERICAL RESULTS 



The input parameters used in the text are given below: tub = 5.279 GcV, ttik — 0.494 GeV |23[. m/, c = 3.524 GcV 
0,S]. For the CKM mixing parameters, we take s 12 = 0.2243, s 23 = 0.04 1 3, s 13 = 0.0037 and S 13 = 1.05 23]. The 
parameters appeared in the wave functions are put in the Appendix B following the corresponding wavefunctions. 

Using the above parameters, finally one obtains the branching ratio of B~ — > h c K~ (in the numerical calculation, 
we take the central values listed in the data book for the input parameters) 

BR{B- -> h c K~) = 3.6 x 10~ 5 . 

Since all the parameters adopted in our numerical computations arc included in the wavefunctions of the concerned 
hadrons and the Sudakov factor which are obtained by fitting data of some well measured processes, thus the error 
of our final result is due to the uncertainties in such fitting. Because of the uncontrollable factors, one cannot expect 
to get very accurate value at the present stage. Therefore, we only keep two significant figures without explicitly 
marking out the error range. But in the last section for discussion, we will present a rough error estimate which may 
make sense. 

There are some theoretical uncertainties in our calculations. One of them comes from the next to leading order 
corrections to the hard amplitudes [24|. In view of this point, we check the sensitivity of the decay rate with different 



choices of hard scales, i.e. we set the hard scales as 

max(0.75^/T c , 0.750k, 1/bi, V&'i, W l/|bi - b q |) < t c < maa:(1.250£ J 1.250k, l/h, l/b[, l/b q , l/|bi - b q |) 
max(0.75 04^, 0.750^, l/6i, 1/6'x, l/b q , l/|bi + b q |) < t d < max(1.25^A^ : 1.250k, 1/bi, l/b[, l/b q , l/|bi + b q |), 

and other parameters are fixed. Then we can obtain the branching ratio of B~ — > h c K~ and the error may be a few 
percents and the whole result is not sensitive to the change of hard scales. 

Another uncertainty comes from the non-perturbative parameters in meson wavefunctions, such as u>b in B meson 
wavefunctions, although they are determined directly from previous experiments or some non-perturbative methods 
like QCD sum rules. Here we vary the value of uib at the range 0.32 ~ 0.48 GeV, with other parameters fixed, as Ref. 



24| did, then we find that the errors can be as large as 20%. It can be observed that the decay rate is relatively more 



dependent on the value of Ub- 
In Ref. Q, 

the author proposed that h c — > r) c j is a promising mode to look for h c and presented a theoretical 
estimation on the branching ratio as BR(h c — > 7? c 7) — 0-50 i 0.11. Combining this branching ratio with BR(j] c — > 
KK-k) — (5.7 ± 1.6)%, one could predict the branching ratio of the decay chain B~ — ► h c K~ — > (r) c ^)K~ — > 
((KKn^K- as 

-> /i c /v- -> (ricl)K~ -> ((KK7r)-f)K-) = 1.0 x 10~ 6 . (7) 
IV. DISCUSSION AND CONCLUSION 

In this work, we employ the framework of PQCD approach to calculate the transition rate of £>~ — > h c K~ and 
predict the branching ratio to be about 3.6 x 10 . 

It is observed that only non-factorizable diagrams contribute to the decay amplitude, because this transition is a 
G-parity violation pro cess. So one can expect it to be relatively suppressed compared with B~ — > J/ipK~ which is 



about 1.00 x 10~ 3 [23]. In our work, we estimate the ratio and obtain results which are reasonably consistent with 
data. In this way, we not only naturally reproduce the production rate of h c at Babar and Belle, but also have a 
chance to study the non-factorizable diagrams. Usually, in most processes, both factorizable and non-factorizablc 
diagrams contribute and it is hard to separate their contributions. Thus when comparing with data, there is an 
obvious uncertainty. However, this case is an optimistic place where only non-factorizable diagrams contribute. It 
enables us to uniquely investigate the non-factorizable contributions to the weak transitions and testify applicability 
of the theory, PQCD. 

Our prediction is consistent with the upper bound reported by the Belle and Babar collaborations ^] at the level 
of order of magnitude. Thus we can conclude that the PQCD is applicable to the process, especially to deal with the 
non-factorizable diagrams. 

In fact, all these phenomenological parameters in the wavefunctions of B,K,h c are priori determined and the 
choice of the factorization scale follows the conventional way. In our numerical computations, we do not introduce 
any free parameter to adjust. Therefore the error of the final result obtained in this work is due to the uncertainties 
included in the wavefunctions of the concerned hadrons and the Sudakov factors which are expressed in terms of a 
few phenomenological parameters. 

As the input parameters which exist in the wavefunctions of the concerned hadrons, vary within a reasonable range, 
an error of about 20% is expected. On other aspect, it would be difficult to make a precise estimate on the uncertainty 



because it is determined by non-perturbative QCD effects, thus only the order of magnitude of the results can be 
trusted. Even then, our theoretical predictions are more consistent with the present data, i.e. the upper bounds set 
by the experiments. It seems that based on the PQCD approach, the obtained results arc more consistent with data 



than the earlier analyses 



12 



3, 



16j, thus one expects that the theoretical framework may reflect the real physics 
picture, at least works in this case. 

If the future experiments further reduce the upper bound and the theoretical estimate on the branching ratio of 
B~ — > h c + K~ cannot tolerate the new data even considering the range for the input parameters, one would confront 
a serious challenge to our understanding of the physics mechanisms, therefore, we hope that by improvement to reduce 
both statistical and systematic errors, new measurements of the Babar and Belle collaborations will provide important 
information towards the answer, otherwise we may need to wait for the data of the future LHC-b. 

To be more accurate in both theory and experiment, we need to wait for more data and improvement in the method. 
Indeed, a favorable channel for observing h c is decay of i(i(2S). However, unfortunately, the BES collaboration has 
not observed this mode and one needs to wait for the upgraded stage, i.c BES III which would accumulate much 
larger database in future and then we can further investigate the consistency between theory and experiments. 
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Appendix A: Some relevant functions appeared in the text 

(i) 

The explicit expressions of decay amplitude M\ appeared in the text can be written as 



M\ 



(c) 



dxx J dx[ I dy I hdh J b q db q I d9[V cb V* s C 2 ((i) - VaVZidbi) + -e c CM)]^ B ^iM)4>l a {v) 

a s (t c ) exp[-S(t c )]n c (x ll x' 1 ,y, b±,b q ), (8) 



(2iV c )3/2 



M. 



(c) 



-16 



dxi I dx[ I dy j bidh j b q db q j dff[C 6 (fi) + -e c Cg(n)](-V tb V t * s )(/)B(xi, h)(f)\ c {y) 



(2iV c ) 3 / 2 



a s {t c ) ex$[-S(t c )]U c (x\,x! x ,y,b\,bq), 



(9) 



8 



M[ d) 



(2Ay 3 / 2 
4Wi 



d Xl fdx[Jdy Jbxdbx [bgdbgJdBWM&M-ViblftiCiM 



1 



(</) 



a s (t d ) e^p[-S(td)]0,d(xi,x[,y, h,b q ), 



(10) 



M. 



-16 



^ 2N ^3/2 J dx x J dx[ J dy j b x db x J b q db q J d6[C 6 (fj.) + -e c C s (fj,)](-V tb V t * s )(t)B(xi, h)(/>l e (y) 



a s (td) exp[—S(td)]£ld(xi, x[, y, bi,b q ), 



(11) 



where the explicit expressions of KA" , KF~ and KT t which come from the contraction of hard kernel and hadronic 
wave functions are given below. Here the meaning of indices i, j have been shown in the text. 

(12) 
(13) 
(14) 
(15) 
(16) 
(17) 
(18) 
(19) 
(20) 
(21) 
(22) 

The explicit forms of £l c ,d(xi,x'i,y,bi,bq) which come from Fourier transformation to products of propagators 
corresponding to quark and gluon are listed as follow 



KA\ c) 


= -16m&V(-l + r 2 )(-.Ti + x[ - 


y), 


kf[ c) 


= —16mb 3 mQr 2 (—xi + x\ — y), 




KT[ c) 


= 16mb 3 m^r 2 (—xi + x\ — y), 




KA 2 c) 


= KA[ d) = 0, 




KP 2 c) 


= —8mb 3 m r 2 (—xi + x\ — y), 




KT 2 c) 


= -8m6 3 m^r 2 (-xi + x\ - y), 




KP^ 


= 16mb 3 rriQr 2 (-l-x 1 +x' 1 +y), 




KT[ d) 


= 16mb 3 niQr 2 (-l - x x + x[ + y), 




KA 2 d) 


= 8m6 4 r 2 (-l + r 2 )(-l-a;i+a;' 1 - 


\-v), 


KP 2 d) 


= 8mb 3 mQ (— 1 — X\ + x\ + y), 




KT 2 d) 


= -8mb 3 mo(-l-x 1 +x' 1 +y). 





n c (x 1 ,x' 1 ,y,b 1 ,b q ) 



n c (x 1 ,x' 1 ,y,b 1 ,b q ) 



[K a {^A c \h q \)e{A c ) + |[-iVo(VXim) + iJ (VA~ c \b q \)]e(-A c )} 

xi^o(V^lbi-M), 

{K {^fAd\h q \)e{Ad) + \[-N Q (^A~ d \h q \) + U {^fAd\h q \)]6{-Ad)} 
xiv ( v ^|b 1 +b,|), 



(23) 



(24) 



9 



with 



A c = m 2 c - Oi + y - 1)((1 - x[)(l - r 2 ) + yr 2 )m 2 b , B c = (1 - Xl )(l - x[)(l - r 2 )m 



A d = m 2 c ~ ( Xl - y)((l - x[)(l - r 2 ) + (1 - y)r 2 )m 2 , B d = (1 - Xl )(l - x[)(l - r 2 )m 2 . 

Here Ji, iVj are 'ith' order Bessel functions of first and second kind respectively, Ki denotes 'ith' order modified Bessel 
functions. 

The explicit expressions of Sudakov factor coming from the resummation of double logarithm appeared in high 
order radiative corrections to the diagrams are also given below 



S(t c ) 



S(t d ) 



- x 1 )p+,b 1 ) + s(x' lP '-, b[) + s((l - x'Jp'-X) 

-\n{t c /A QCD ) t ^ -ln^/Agoo) , ^ - ln(i c /A QC £>) 



1 

Ti 



In. 



In. 



-\n{biK QC D) Iu^Aqcd) \n(b q A QC D) 

s((l - a;i)p + ,6i) + s(x[p'~ ,b[) + s((l - x'^p 1 ' ,b' x ) 

-]n{t d /A Q c D ) , ^-Ht d /A QCD ) t - hx(t d /A QC D) 



1 



In 



ln(6iA QC r>) ' ^ ln(b[A Q cD 



ln{b q A QC D) 



(25) 



(26) 



where 



t c = max{^/A~ Cl y/B~ c , 1/h, l/b' v l/b q , l/|b x - b q |), t d = max(y^, y^B~ d , 1/h, l/b[, l/b q , l/|b x + b q |). 



The explicit expressions for the Sudakov factors are given in Rcf. 



26J as 



s{uj,Q) 



A = C F ^ 



hx(Q)A[a a (p)]+B[a a (p)] 



7 9 (a«(/i)) = -a a ((j,)/w, 
_ 33-2n f 



"67 7T 2 10 8 „ , ,eT B , 

nf H — Pnlni 

9 3 27 } 3^ v 2 ; 

=27e-1 



(-) 2 , 
7T 



where 7b is the Euler constant, nf is the flavor number, and 7 g is the anomalous dimension. We will take nf equal 
to 4 in our numerical calculations. 

Appendix B: Wave functions relevant to calculation 



The B meson light cone wavefunction is usually written as 



27] 



1 ^4 



(27T) 



Z e* k - z (0\q p (0)b a (z)\B(p)) = -^={(i, + mj3 ) 75 ^ B (k) - ^— ^ B (k) 



'2N C 



V2 



(27) 



a/3 



10 



where n = (1,0, Ot) and v = (0, 1,Ot) denote the unit vectors corresponding to the "plus" and "minus" directions 
respectively. In eq. (|2T|) . two different Lorentz structures exist in the the B meson wavefunctions. <pB (k) and <pB (k) 
satisfy the following normalization conditions respectively 



0s (k) 



h 



d 4 k 



{2ir)^ ay ' 2V2N~ C ' J (2tt) 4 

In the numerical calculation, one usually ignores the contribution of </>s(k) 
from 



s(k) = 0. (28) 
and only take the contribution 



28 



1 



/2N, 



:(^ + ms)750s(k). 



(29) 



The wave function of B meson is given as 

4>B{x,b) 



28 



--fBX 2 (l - x) 2 exp 



2V27V C 

where uo b = 0.4 GeV and N B = 2.4 x 10 3 . The decay constant of B meson f B = 0.18 GeV. 
The twist-3 light cone distribution amplitude of K meson is expressed as 

fi 



(30) 



(K°(j/)\a a {z)qp(0)\0) = 



/2N, 



where ruff 



given as 



c J o 



dxe^ -{ 75 j/<£(x) + mX^KW + <Mjb A - WkW} 0a , (31) 



m s +mi 



In a recent work 



30] . the K meson wave function distribution amplitudes used in eq. (|3 1 1) is 



2V2lV^ 



:{6s(l-x) (l + afC? /3 (2a:- l) + a%Cl ,2 {2x- 1))}, 



(32) 



: {1 + 3p£(l + 6of ) - 9p A 'af + C 1 1/2 (2a; - 1) 



2J2N C 



2' if if K ( % „_ A' 

yP+ a i -P- I g + 27a 2 



+C 2 1/2 (2x - 1) (30 mK + 15^a? - 3 P K a? ) + C 3 1/2 (2x - 1) fl0 % if A 3A - - | p5af 



-3r) 3K uj 3K Cl /2 {2x-l) + ^(p$ + p*)(l-3af + 6af)ln: 

+ |(p^-/^)(l + 3of + 6af)lnx}, 

iT _ 1 #M ^ 
A 6 da; ' 
/if 



2x/2lV^ 



{6a;x 



1 + § A? + 15p^af - H p*af + Up* a? y « ) C?' 2 ^ 



+ (5r?3if - ^ V3KW3K + \ Cl /2 {2x - 1) + V3K X 3K C 3 3 /2 (2x - 1) 

+ 9xx(p*f + p A ')(l - 3af + 6af ) Ins + 9;rx(P+ - /9 A ')(1 + 3af + 6af ) lnx} . 



(33) 
(34) 



(35) 
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with 



K _ (m s + m q ) 2 

P+ ~ 2 : 



K _ m 2 s - m\ 



K 



V3K 



Jzk m q + m s 



Ik m 2 K 



Cl /2 (t) = t, C* 2 1/2 (i)^(3i 2 -l), C 3 1/2 (t) = i(5t 3 -i), C 4 1/2 (i) = i(3-30t 2 + 35i 4 ), 
Cl'\t) = 3t, C 3 2 /2 (t) = ^(5i 2 - 1), Cl' 2 = i(35i 3 + 3t), 



where af = 0.06 ± 0.03, 



0.25 ± 0.15, / 3 k = (0.45 ± 0.15) x 10~ 2 GcV 2 , lo 3 k = -1.2 ± 0.7, A 3 *r = 1.6 ± 0.4, 



/k = 0.16 GeV, m s = 137 ± 27 MeV, m q = 5.6 ± 1.6 MeV. 

The light cone distribution amplitude of h c meson is proposed in 



20] 



(^(5,6)1^(^(0)10) = 



<2N, 



due 1 



c JO 



m\ c fee* ■ z 



f ho m hc 4\\ + , , 2 75<(") - ft Al Al^l(u) (36) 

L Z \1 ' Z ) J J /3a 



with 



q-z 



9n 



q-z 



We have dropped out the twist-3 distribution amplitudes here, because h c is heavy, the higher twist contributions 
are negligible. Because the produced h c in the transition B~ — > h c K~ can only be longitudinally polarized, the 
wavefunction <j>j^ (u) does not contribute and we neglect its explicit form in the text. It is a good approximation to 
assume that (/)» (u) is of the same form as Xo which is the leading- twist distribution amplitude of \c0 defined in Eq.(59) 



of Ref. 



19] 



4>{ (x) = 27.46- 



fhc 



:(1-2X) 



x(l - x)[l - 4x(l - x) 



2^J2N C *- '[ [1 - 2.8x(l - x)] 2 

where h c decay constant fh e is set to be 0.335 GeV which is the same as f Xal as has been done in Ref. 



(37) 
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